Directed polymers with tilted columnar disorder and Burgers-like turbulence 
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The minimal energy variations of a directed polymer with tilted columnar disorder in two dimensions 
are shown numerically to obey a multiscaling at short distances which crosses over to global simple 
scaling at large distances. The scenario is analogous to that of structure functions in bifractal 
Burgers' turbulence. Some scaling properties are predicted from extreme value statistics. The 
multiscaling disappears for zero tilt. 
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The directed polymer in a random medium (DPRM) is 
a prototype model in the statistical physics of disordered 
systems [Q. At zero temperature, T, for each realisation 
of a random potential landscape, the problem consists in 
the statistical characterisation of long directed paths (no 
overhangs) which minimise the total energy. Besides hav- 
ing ramifications and generalizations in different fields, 
the DPRM is directly related to flux lines in type II su- 
perconductors, domain walls in random fcrromagnets and 
nonequilibrium fluctuations of growing interfaces in di- 
mension d = 2 ||. In the last case the DPRM energies 
can be mapped onto the interfacial profiles realised dur- 
ing the growth, and the quenched average over potential 
disorder amounts to a mean over growth histories. 

Attention has been devoted in recent years to the 
effects of various forms of disorder on the asymptotic 
DPRM scaling properties. Particular interest focused on 
long-range disorder correlations which in the language of 
interfacial growth, can be cither in "time" t (longitudi- 
nal spatial variable) or in space x (transversal) S, |Q. 
While the effects of spatial correlations are relatively well 
understood, long range temporal correlations (including 
those appropriate to 1// noise) remain more problem- 
atic mainly due to a lack of Galilean invariance. Here, 
we study the case of a noise which is infinitely correlated 
along one space-time direction (x — — t), a situation we 
refer to as tilted columnar disorder In the DPRM 

context, columnar disorder is of particular interest for de- 
scribing flux lines in superconductors irradiated by heavy 
ions Q. The case in which the columns are tilted with 
respect to the (time) direction of the flux lines could be 
of relevance, e.g., for the behaviour of a flux line in the 
presence of a forest of splayed columnar defects. Such 
situations have been studied recently both theoretically 
and experimentally ^j^]. 

In this Letter we show that the DPRM with tilted 
columnar disorder in d — 2 and at T — is characterised 
by two distinct critical regimes. The first is a global one, 
holding at large length scales and consistent with sim- 
ple scaling. We fully elucidate this regime in terms of 
extreme value statistics. The second regime is local and 



FIG. 1. Tilted columnar defects are oriented along lines 
x + t constant. At T = 0, the optimal paths (thin lines) follow 
columns with very low energies £q,£i,£2,-- (thick lines) for 
a very large fraction of time (see text). Periodic boundary 
conditions are applied at the horizontal sides of the strip. 



holds at smaller length scales. It is characterised by in- 
termittency and multifractal scaling. The two regimes 
match in a full scaling framework which displays close 
analogies and correspondences with the scaling regimes 
of a fully turbulent flow p|-|lC|] . The very existence of this 
unexpected local multiscaling is especially remarkable be- 
cause the distributions generating it are quite standard 
and general, and do not necessarily involve long tails. 

Consider a directed polymer that evolves on a 1 + 1 
dimensional lattice in a random energy landscape e(x, t). 
The energies e(x, t) associated to each site (x, t) are taken 
from a probability distribution function (pdf) P{e). The 
partition function Z(x,t) is then defined through the re- 
cursion relation 

Z(x, t) = e~ e{x > t)/T [Z{x - 1, t - 1) + Z(x + 1, t - 1)] (1) 

together with an initial condition (we take Z[x, t = 0) = 
e -e{x,o)/T w jt n £ [x^ o) randomly distributed according to 
P(s)). We will in particular study tilted columnar de- 
fects, e(x,t) = e(x + t), and will compare them to the 
more standard columnar defects, e(x, t) — e(x) j^j. Fig. 
1 sketches the situation for the tilted case. 

It is generally assumed that the recursion (|l]) is a 
discrete version of the diffusion equation 



dZ_ 
~dt 



(2) 
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where V(x, t) is a random potential with zero mean. 
From (||) one obtains the KPZ equation for the height 
h of a growing interface by the usual Hopf-Cole transfor- 
mation h = T log Z 



dh _ Td 2 h 
dt ~ ~2~£h? 



dh 

dx 



(3) 



Finally, (y) can in turn be mapped onto the Burgers equa- 
tion Q 
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9ii 
9t 



<9w 
eta 



(4) 



by letting u = —dh/dx, where K is interpreted as the 
velocity of a fluid subject to a random force rj(x, t) = 
—dV/dx, and with viscosity T/2. 

Here, we will limit ourselves to the T — ► limit of (Q). 
Denoting limy_^ l°g ^ = — -^j (0) reduces to 

J5(x, i) = min[E{x - 1, t - 1), E(x + 1, t - 1)] + e(x, *) (5) 

Taking s = in (|) one obtains, in the continuum limit, 
an equation with a nonlinearity \dh/dx\ Noise dis- 
order generated when e ^ is however expected to trans- 
form such nonlinearity into the analytical one of (|^). Us- 
ing the numerical methods of rehj lj] , we found evidence 
of such analytical nonlinearity p5 |. 

In Ref . || , some properties of the minimal energy path 
in columnar disorder (t? = e(x)) were studied numerically. 
It was found that the disorder averaged transversal dis- 
placement R(t) = (x m [ n (t)} °^ the mm i ma l energy path 
grows as a power of time R(t) ~ t 1 /*, where z is highly 
non-universal, and can be calculated for various pdf's 
P using Flory arguments based on extreme value statis- 
tics. In particular for a Gaussian pdf, z = 1, but with 
logarithmic corrections. These arguments should remain 
valid for the tilted columnar disorder we study. Indeed, 
z remains unchanged as we verified numerically for a few 
distributions P(e). 

Here, we study the moments of the energy-energy 
structure function 



G q (r,t) = (\E(x + r,t)-E(x,t)\i) 1 ^ 



(0) 



and find that for the case of tilted defects, and for a 
very large class of pdf's P{s), G q (r, t) shows two scaling 
regimes. The first one is global for r > R(t) ~ t 1 / 2 , with 
G q (r, t) independent of r and with g-independent simple 
scaling in time 



G q {r,t) ~R(t) a ~t 



a. i a I z 



> R(t) 



(7) 



In contrast, for r < R(t), there is strong evidence of 
multiscaling 



G q (r, t) ~ r c « 1 < r < R(t), t fixed 



(8) 
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FIG. 2. Log-log plot of G q (r,t = 2048) versus r for 
= 1, ... ,7. The data are for P(e) = 5|e|" 6 with e < -1 




FIG. 3. Log-log plot of G q (r = 1, t) versus t for 
q — 1, . . . , 7. The data are for a standard normal distribu- 
tion. 



Figure 2 shows our numerical results for P(e) = 5|e| -6 
with e < — 1. These results were obtained from a nu- 
merical iteration of (||) on a strip with periodic boundary 
conditions and a typical width 2 15 in both space and time 
directions. The results were averaged over 2 x 10 3 realisa- 
tions of the random energy landscape. Qualitatively sim- 
ilar behaviour was found for a Gaussian, and for power 
law pdf's like P+(e) = /ie _M_1 (e > 1) or even for pdf's 
with a bounded support like P v (e) = (y+ l)e" ', e € [0, 1]. 

We also investigated in detail the nearest neighbour 
energy difference which, when replaced by a space deriva- 
tive, should correspond to the local velocity u in (Q). Also 
this quantity shows multiscaling, i.e. 



(\<u\«) 1 l«{t)~G q {r = l,t)~t 



a q /z 



(9) 



where a q is g-dependent. In Figure 3 we show our results 
for for a Gaussian distribution. Qualitatively 

similar behaviour was found for all P's investigated. 

The scaling behaviour of G q (r,t) as given in 
and (|^) can in fact be described by a form which was 
originally proposed for the structure functions of fully 
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Tabic 1: Numerical estimates of exponents for P 2 {s) = 3e 2 
with e £ [0,1]. The error in these exponents is estimated 
to be a few percent. 
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developed turbulence || , and later also for height corre- 
lators of certain models of epitaxial growth fin] : 



G q (r,t)~R(t) a *r^f q 



R(t) 



(10) 



where f q (u) — > constant for u — > 0, and f q (u) — > u~^ q 
when u — > oo. Consistency with (^) then implies the 
scaling identity 



(11) 



In the scahngs (jq)-([ll|) the DPRM minimal energy 
E(x,t) plays the role of a fluid velocity component in 
turbulence, while R corresponds to the Reynolds num- 
ber. Moreover, [G 2 (r = 1:^)1 2 is an analog of the local 
average energy dissipation || [l(J. In all cases investi- 
gated, (|l|) holds within the numerical accuracy. As an 
example, Table 1 shows our estimates for a q , £ q and a for 
some q- values for the case of P 2 {e) = 3e 2 with s G [0, 1]. 
Since in all cases q( q ~ 1 for all q > 1, while roughly 
qC, q ~ q for g < 1, the multiscaling appears consistent 
with a biscaling, with £ q very universal with respect to 
the noise distribution. On the other hand, as discussed 
below, a is less universal, and its classes appear to be at 
least as many as the possible asymptotic distributions in 
extremal value statistics. Biscaling is a form of multiscal- 
ing already detected for velocities in Burgers turbulence 
p6[ . We will give full details of all exponent estimates 
for different P's elsewhere. Eq. ( p0| - |TT| ) reveal that, in 
presence of a tilt, the variations of minimal energy corre- 
sponding to even relatively small r are wildly fluctuating 
in the same way as longitudinal velocity differences in a 
turbulent flow. 

The geometry of minimal energy paths suggests an an- 
alytic determination of a. Consider the minimal energy 
path starting somewhere and arriving in xq at time t (Fig. 
1). Typically, this path will have followed a column with 
a favourable low energy £o for a 'macroscopic' fraction of 
t. The particular path chosen will depend on a balance 
between the energy gained there and the energy it costs 
to get from the starting point to the column of energy 
Eq and then to xq in time t. As time goes on it pays to 
'travel' to more distant but energetically more favourable 
columns. Now consider paths arriving at x\ ^ Xq. For 
small \xq — these paths will in general have followed 
for a large fraction of t the same low energy column. 



Thus, the difference in energy between paths ending at 
xq and at x\ will arise in a short part of their histories, 
where they did not follow the same column. Considering 
sites X2 that are further and further away from ieo, we 
will arrive at the situation where it is more favourable 
for the path ending at X2 to follow another low energy 
column with energy e\, closer to x 2 , for most of its his- 
tory (Fig. 1). The typical distance between the different 
minimal energy columns must be of the order of R(t), 
the only relevant length scale in the problem. For the 
most part, each path follows the column with the lowest 
energy within a local region of size R(t). Subsequently, 
at scales above R(t) the energies of different paths are 
essentially independent, which accounts for the constant 
behaviour of G q on these scales visible in Fig. 2. As 
explained above, the columns that attract most of the 
paths in a region of size R are such that their energy is 
minimal in a region of that size. The minimum e in a set 
of R independent random variables drawn from P(e) is 
distributed according to 



P R (e) = RP(e) 



1 - 



P{e')de' 



R-l 



(12) 



Since for r > R(t) energy differences between two paths 
are, up to small corrections, given by the energy differ- 
ences of the minimal energy columns they follow, and 
since these are expected to be independent, we may write 

[G,(r,t)]« J J \e 1 ~e 2 \ q P R (e 1 )P R (e 2 )de 1 de 2 (13) 

It is not possible to evaluate (|l^) for arbitrary pdf 's an- 
alytically. However, from extreme values statistics [jl7j it 
is known that for R — > oo, P R converges to a limited num- 
ber of possible asymptotic forms which only depend on 
the behaviour of P(e) at — oo. As an example, for a Gaus- 
sian P, or for other pdf 's that approach zero at — oo faster 
then any power, Pr(s) converges to the Gumbel pdf 
Pg{v) = exp(y-expy) where y = V2 log R(e+ y/2 log R) . 
Inserting this into (p~3) one finds 



[G 9 (r,i)r (log i?)-«/ 2 



(14) 



Since for the Gaussian distribution R(t) ~ t/logt [EJ, 
G q (r, t) indeed satisfies (0), with a — 1, up to logarithmic 
corrections. One can proceed similarly for other pdf's. 
For distributions whose support is bounded from below 
such as P + , or P v (e), Pr converges for R — > oo to the 
Weibull distribution, from which the universal value a = 
1 is again obtained, but without logarithmic corrections. 
Finally, for pdf's approaching zero as a power law at — oo 
, i.e. P-(s) ~ a is non universal and equals 1 + 

l//x. These predictions are consistent with our estimates 
a - 1.03 (P+(e) with n = 3), 1.05 (P 2 (e)) and 1.19 ( 
P_(e) with fj, = 5). For the Gaussian case the data can 
be fitted well by the form (Q). 
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We verified that when the tilt is removed and when 
e(x, t) — e(x) is assumed, the multiscaling disappears for 
all P's studied. However, a remains unchanged. Indeed, 
our arguments for the its determination should be valid 
both in presence and in absence of tilt. 

The difference between these two cases can be under- 
stood by referring again to Fig.l. Let us assume that e\ 
is particularly low. Only paths whose endpoint is above 
that column can take advantage of this low energy. As 
a consequence, there is a very big energy difference for 
optimal paths ending at sites x that are close to, but on 
different sides of the column with energy e\. In contrast, 
if the tilt is absent, two optimal paths ending on different 
sides of a very low energy column can both follow that 
column for a long time, and hence energy differences tend 
to be smaller. The large energy differences in the tilted 
case lead to broad distributions and multiscaling. 

It remains an interesting question to find out whether 
the multiscaling observed here is present for any value of 
the tilt angle, or whether it only appears above a certain 
non-zero critical angle. 

In ||, the KPZ equation (||) with noise V(x,t) that 
has long range correlations in time 

(V(x, t)V(x', t')) ~ \t - t'\ 2e - l 8(x - x') (15) 

was studied. From a renomalisation group analysis it was 
expected that for 6 sufficiently below 1/2 

(l + 20> = 2a. (16) 

The case 8=1/2 marks the transition to a regime where 
infrared divergences become important and the renor- 
malisation procedure breaks down. Dimcnsionally, the 
noise corresponding to columnar disorder in absence of 
tilt, shares the same correlation ( |l5| ) with the l/f noise 
(9 = 1/2). Thus, once accepted that (g) at T = is the 
appropriate continuum limit, our results suggest that for 
pdf's (e.g. Gaussian) whose extreme value statistics is 
given by the Gumbel distribution, ( |l6| ) remains valid at 
8 = 1/2, but with logarithmic corrections. 

In summary, the differences in optimal energy of the 
DPRM with tilted columnar defects obey a form of mul- 
tiscaling analogous to that valid for velocity structure 
functions in turbulence. An intermittent scaling regime 
matches a global simple scaling one. The global regime 
can be characterized on the basis of extreme value statis- 
tics arguments. The multiscaling at small scales is an 
unexpected feature of both deterministic and stochastic 
models of nonlinear diffusion. The multiscaling of the 
discrete counterpart of the velocity potential observed 
here is similar to the scaling behaviour of velocity incre- 
ments of the Burgers equation solutions with other types 
of random forcing [|l6| . The case treated here is new also 
because, if (^) indeed provides the appropriate contin- 
uum limit, the biscaling is valid for increments of the 
velocity potential. 



The research presented here can be extended in several 
directions. Besides obvious generalisations to finite tem- 
perature and higher dimension, the case of networks of 
splayed columnar defects deserves special attention [pj. 
It would be worth investigating whether the multiscal- 
ing observed here also shows up in these experimentally 
relevant situations. 
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